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Modeling a single crystal of cuprate high- Tc superconductor, such as Bi2Sr2CaCu208+4, as a stack 
of intrinsic Josephson junctions, we formulate explicitly the cavity phenomenon of plasma oscilla- 
tions and electromagnetic (EM) waves in mesas of cylindrical and annular shapes. When the mesa 
thickness is small compared with the EM wavelength, the boundary condition for the inductively 
coupled sine-Gordon equations is the Neumann-type one to a good approximation, addressed first 
theoretically and verified in a recent experiment. This renders the superconductor mesa a cavity. 
Biasing a dc voltage in the c direction, a state with ±7r kinks in the superconductivity phase dif- 
ference piled up alternatively along the c axis is stabilized. The ±7r phase kinks provide inter-lock 
between superconductivity phases in adjacent junctions, taking the advantage of huge inductive 
couplings inherent in the cuprate superconductors, which establishes the coherence across the whole 
system of more than ~ 600 junctions. They also permit a strong coupling between the lateral cav- 
ity mode and the transverse Josephson plasma, enhance the plasma oscillation significantly at the 
cavity modes which radiates EM waves in the terahertz band when the lateral size of mesa is set to 
tens of micrometers. In order to overcome the heating effect, we propose to use annular geometry. 
The dependence of frequency on the aspect ratio is analyzed, which reveals that the shape tailor is 
quite promising for improving the present technique of terahertz excitation. The annular geometry 
may be developed as a waveguide resonator, mimic the fiber lasers for visible lights. 

PACS numbers: 74.50.-|-r, 74.25.Gz, 85.25.Cp 
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I. INTRODUCTION 

The Josephson effect provides a unique principle to ex- 
cite high-frequency electromagnetic (EM) waves [J, 
Much effort has been made to stimulate powerful ra- 
diations, first using artificial Josephson junction arrays 
d, 0, H, ll, 0, B H, and later on Josephson junctions 
inherent in cuprate high- Tc superconductors of layered 
structures [HI, lil, [H, M E H [B M, H S El, [H, 
m, [H, [2^. The latter have obvious advantages, since 
the junctions are homogeneous at the atomic scale guar- 
anteed by the high quality of single crystals, and the 
superconductivity gap is large, typically of tens of meV, 
which in principle permits the frequency cover the whole 
range of the terahertz (THz) band, a very useful regime 
of EM waves still lacking of compact solid-state genera- 
tors Hi, [23. 

An experimental breakthrough was achieved in 2007 
[2^ . Clear evidences have been obtained that coherent 
terahertz radiations from side edges of a thin rectangu- 
lar mesa were realized by biasing a dc voltage in the c 
axis of the Bi2Sr2CaCu208+5 (BSCCO) single crystal; 
the frequency of the EM wave and the voltage where the 
radiation occurs obey the ac Josephson relation, and the 
frequency coincides with one of the cavity modes deter- 
mined by the lateral size of mesa [2^ . 

This discovery is expected to leave significant and long- 
standing impacts. It has the potential to open a practical 
way to develop a source of frequency tunable EM waves 
based on superconductivity and fill the so-called THz- 
gap. It also shines light on new directions of making use 
of the phase of superconductivity. Its importance would 



be better appreciated if one notices that up to now the 
usage of superconductivity phase variable is still limited 
to the SQUID. 

The discovery immediately raises many interesting 
questions, such as why the Josephson plasma oscillation, 
namely the coherent tunneling of Cooper pairs forth and 
back between adjacent CuO layers driven by a c-axis volt- 
age can radiate strong transverse radiations, in absence 
of Josephson vortices induced by an applied magnetic 
field; how it becomes possible to synchronize the super- 
conductivity phase variables of ^ 600 junctions, and so 
on. These quests challenge our knowledge on supercon- 
ductivity as well as nonlinear phenomena. 

It is formulated explicitly in Ref. that there is a sig- 
nificant mismatch in impedance |^out|/|^in| ~ Aem/^z 
at the edge of a superconductor mesa when the thickness 
of the mesa is small compared with the EM wave- 
lenght Aemj which was known for a single junction as 
a limiting case [2^. This makes the tangential compo- 
nent of oscillating magnetic field at the mesa edge ex- 
tremely small compared with the electric one, in sharp 
contrast to the case of EM plane waves. Actually, this 
relation is a general property shared by a normal ca- 
pacitor with small electrode separation. Particularly for 
thin mesas of superconductor, the vanishingly small tan- 
gential component of magnetic field gives the Neumann- 
type boundary condition of superconductivity phase dif- 
ference across junctions, namely the spatial derivative of 
the phase difference normal to the edge should be zero 
(in absence of an applied magnetic field). 

With this Neumann boundary condition, a new dy- 
namic state of superconductivity phase difference has 
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been found in a stack of Josephson junctions under a 
dc voltage bias, which is characterized by static ±7r (ac- 
tually ±(2m + l)7r with m integer) phase kinks localized 
at the mesa center and piled up alternatively along the 
c axis [30I. [31I [3^. [33} . In this tt kink state, an inter-lock 
between superconductivity phases in neighboring junc- 
tions appears, which establishes the coherence across the 
whole system of more than ~ 600 junctions. The tt phase 
kink makes the lateral cavity modes and the transverse 
Josephson plasma couple strongly to each other, permits 
large dc supercurrent flow into the junctions and gener- 
ates strong EM radiations from the sides of mesa. 



The TT kink state requests strong inductive couplings 
between adjacent superconducting layers [13], which 
is guaranteed in the cuprates high-T^ superconductor 
BSCCO. Another condition is on the thickness of the su- 
perconductor mesa: it should be much smaller than the 
EM wavelength on one hand, which renders the supercon- 
ductor mesa a cavity, and on the other hand, it should be 
thick enough to avoid large surface effects. The samples 
of Lz — 1/im adopted in the experiment [2^ satisfy these 
requests well. The tt kink state is stable against other 
distortions, such as thermal fluctuations, inhomogeneity 
of physical parameters, and so on. 



The present work is motivated by two recent exper- 
iments on THz radiations from thin mesas of BSCCO 
single crystals. The first one is a detection of THz radi- 
ations from a cylindrical mesa [33|, after the theoretical 
proposal [31] for testifying the boundary condition of the 
system. The second one is an observation on EM stand- 
ing waves in a rectangular mesa (ssj . an approach being 
able to provide direct evidence for cavity resonance. 



The remaining part of the present paper is organized 
as follows. In Sec. II, we discuss the right boundary con- 
dition for thin mesas from a theoretical point of view. 
Then we point out that the result of the recent exper- 
iment using a cylindrical mesa is to be understood in 
accordance with this boundary condition. In Sec. Ill, 
the radius dependence of the radiation frequency for first 
several modes is presented, accompanied with the spatial 
distributions of the n kink in superconductivity phase dif- 
ference, supercurrent, as well as the EM standing wave. 
Section IV is devoted to a new proposal of using annular 
mesa, which can reduce the Joule heating and enhance 
the leakage of heat as compared with the cylindrical one. 
By checking the frequency variation upon tuning aspect 
ratio, it is displayed that tailoring the mesa geometry 
may be a prosperous way to improve the present tech- 
nique for THz stimulation. Finally summary and per- 
spectives are given in Sec. V. 



II. BOUNDARY CONDITION AND 
CYLINDRICAL MESA 

A. Basic equation and general solution 

The inductively coupled sine-Gordon equations for the 
gauge-invariant phase differences in a stack of Josephson 
junctions with a dc current bias and dissipations are given 
in the dimensionless form as [13, HH [111 

APi = (1 - CA(2))(sinPz + f3dtP, + dfPi - Jcxt), (1) 

under the Neumann boundary condition in the later di- 
rections dnPi = 0; C = A^ft/(s + Z))2 ~ 10^ is the induc- 
tive coupling and /3 = AwacXc/c^/e = 0.02 is the c-axis 
conductivity; A the Laplace operator in lateral direc- 
tions, and A'-^^Qi — Qi+i + Qi-i — '2Qi the second dif- 
ference operator along the c axis. The lateral space is 
scaled by Ac and time by c/Ac-y/e. For more details of 
definitions see Refs. (30l[3ll[33j. 

The experimental observation of the cavity relation of 
radiation frequency [23,[13] indicates that standing waves 
of plasma oscillation has been built in the cavity formed 
by the mesa, which in turn implies that the oscillating 
part of the pha se difference satisfies the Laplace equa- 
tion [13, [31I l33l| . The plasma oscillation should be uni- 
form since the observed radiations are coherent, known 
as a superradiation [2^. Without losing generality, the 
solution to Eq. ^ can be given by 

Pi{r,t)=ujt + Pir,t) + fiP'{r), (2) 

where the first term accounts for the finite dc bias volt- 
age, and the second for plasma oscillation 

P{r,t) = Ag{r)siii{ujt + ip) (3) 

with A the amplitude, g{r) an eigenfunction of the 
Laplace equation with the Neumann boundary condition 
for the corresponding geometry, and the frequency given 
by the voltage following the ac Josephson relation; the 
third term carries the inter-junction coupling via the I 
dependence. 

The general form of Eq. ^ describes a wealth of solu- 
tions even giving the constraint imposed by available ex- 
perimental results. However, it is easy to see [13, [Ml [13| 
that the sequences fi = (—1)' or fi = (— with pe- 
riod 2 or 4 layers respectively diagonalize Eq. ^ due to 
the property of the difference operator A*^^^. In simula- 
tions, these two states appear frequently and are stable 
[13, [l^. Here we focus the attention to these two cases. 

The static phase P*(r) should then satisfy the differ- 
ential equation 

AP" = qC cos (pJi{Ag{r)) sin P", (4) 
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with the same boundary condition of the total phase dif- 
ference, where g = 4 and 2 for the c-axis sequence of 
period 2 and 4 respectively. 

The static phase factor P''(r), amplitude A, phase shift 
in plasma oscillation ip and total current Joxt can be 
solved for a given voltage (or equivalently w) [s^, l3l| . 
The nontrivial solution P*(r) presumes the tt phase kink 
[30, 31, 32, 33]. The static c-axis supercurrent which con- 
tributes to the net current flowing into the system in the 
c direction is given by 

J,{r) = sin ipJi{Ag{r)) cos(F"(r)). (5) 

We point out that the claim in Ref. [37j on a very non- 
linear IV characteristics in the absence of tt kinks is not 
compatible with Eq. ^ and needs further check. 

The trivial vacua solution P" = is dropped here, 
which gives a linear IV and a broad voltage regime with 
weak radiations, inconsistent with experiments psi . Isg! ] . 
See Ref. [H, for more discussions on this vacua state. 

It is noticed that only the fundamental frequency is 
taken in Eq. Q since higher harmonics are small es- 
pecially in cylindrical and annular geometry, where the 
higher harmonic frequencies appearing in the expansion 
of the sine function of the dc Josephson relation are not 
the cavity frequency. 

The oscillating electric and magnetic fields inside the 
junctions are given by the plasma oscillation term of the 
phase difference in dimensionless forms as 



E^{r,t) = —, (6) 

with E'^jEy ~ to a good approximation in cuprates 
[lOl, and 

B{r,t)^-Vx{P{r,t)z), (7) 

with z the unit vector in the c axis. It is Eq. ([7]) that 
enforces the boundary condition for the phase difference 
d,iPi = with n the normal of the sample edges. 



difference inside the junctions [sl*]. For the cylinder ge- 
ometry, the radial part of an eigenfunction of the Laplace 
equation is given by the Bessel function, not their deriva- 
tives. The boundary condition for the cylinder geometry 
determines uniquely the wave number, and vice versa, 
since the zeros of the Bessel functions are different from 
the zeros of their derivatives, contrasting to sine and 
cosine functions in the rectangle geometry. Therefore, 
measuring by experiments the frequency of EM radiation 
from a cylindrical mesa of given radius and assigning the 
mode by the properties of Bessel functions enable one 
to identify uniquely the right boundary condition for the 
EM waves. It was addressed that, from the boundary 
condition suitable for thin mesas, the frequency of radia- 
tion should be given by the zeros of derivatives of Bessel 
functions [3ij. A case of the (0,1) mode, which is uniform 
azimuthally, has been worked out explicitly, and the fre- 
quency is figured out as — Xoi/'^ with the radius a 
and Xoi — 3.8317 the first zero of derivative of the Bessel 
function Jo{z) j39| . 

The recent experiment [3^ showed that the frequency 
of EM wave radiated from the cylindrical mesa of radius 
a ~ 45//m is / ~ 0.474TIIz. With the fight velocity in the 
sample c' = c/^e ^ 7.16 x lO^m/s (e ~ 17.54 [ifl]), the 
wave number observed in the experiment [3^ is given by 
k ~ 1.87/a. Since x = 1.8412 is the first zero of deriva- 
tive of the Bessel function Ji{z), it is inferred that the 
mode (1,1) has been achieved in the experiment. The 
most important message from the experiment is that the 
frequency of EM radiation is determined by the zero of 
derivative of the Bessel function, in agreement with the 
theoretical analysis |31]. The right boundary condition 
for the cavity therefore is that the tangential component 
of the oscillating magnetic field should be zero at the 
edges as discussed above, consistent with the analysis in 
Ref. [l^l- The (1,1) cavity mode for the cylinder geome- 
try assumes a lower frequency than the (0,1) mode due 
to the positions of zeros of derivatives of Bessel functions 
[39| . and thus is realized at a low voltage while the (0,1) 
mode might be difficult to excite in experiment where the 
heating effect is crucial. 



B. Testifying the boundary condition 



In a rectangular mesa, the spatial part of the plasma 
term of the lowest cavity mode satisfying the right bound- 
ary condition is cos{xtt/L). Unfortunately, cos{xtt/L) 
and sin(a;7r/L) are both eigenfunctions of the Laplace 
equation, derivative to each other, and give the same 
wave number (or equivalently frequency in dimensionless 
form) tt/L. Therefore, from a pure experimental point 
of view, the dependence of frequency on the system size 
cannot determine uniquely the mode. 

A mesa of cylindrical shape was proposed to testify 
uniquely the dynamics of the superconductivity phase 



III. MODES FOR CYLINDRICAL MESA 



Knowing the right boundary condition, it is not hard 
to figure out the EM standing waves for all the cavity 
modes. Since their distribution can be observed directly 
in experiments [sHi and serves as a check of the theory, 
we map them out explicitly. The distributions of super- 
conductivity phase difference and the supercurrent are 
helpful for understanding the way how large dc powers 
are converted to high-frequency transverse radiations. 
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FIG. 1: (color online). Spatial distribution of (a) the static 
phase term P" , (b) supercurrent, and the standing (c) elec- 
tric and (d) magnetic wave for the (1,1) mode of a cylindrical 
mesa. Here = qA^sinip — 5000 is taken for Eq. Q in (a) 
and the kink is approximated as a step function in (b) corre- 
sponding to = oo. The lateral coordinates are normalized 
by the radius of the cylinder. The quantities except the phase 
difference P' axe up to the plasma amplitude A. 

A. (1,1) mode for cylindrical mesa 

All the modes (m, n) for the cylindrical geometry can 
be classified according to the number of symmetry axes 
in the azimuthal direction, m, and the number of maxima 
and minima along the radial direction, n, in the spatial 
part of plasma oscillation and electric field (see Eq. ([6])), 
which all correspond to the zeroes of the derivative of 
Bessel functions. The solution to Eq. ([T]) assigned as the 
(1,1) mode for the cylinder geometry is given by 

Pi{r,t) ^Ljt + AJi(^^p^ cos(t)sm{ujt + ip) + fiP'{r), 

(8) 

where the cylindrical coordinates are given by r (p, (jj) ; 
Xii — 1.8412 is the first zero of the derivative of Ji{z). 
For a cylinder with uniform physical properties, the 
eigenfunction for the azimuthal angle sin^ degenerates 
with cos(/), which, and any linear combination of them, 
can be absorbed into the latter by redefining the az- 
imuthal angle </>, and will not be discussed explicitly. 

The wave number, and equivalently frequency (k = f 
in dimensionless form and k = c' f with units), of this 
mode is given by k^^ = Xii/a — 1.8412/a, which agrees 
well with the experiment [3J| as discussed in the previous 
section. 

The distribution of the static phase P*, supercurrent 
and the spatial part of the oscillating electric and mag- 
netic fields are displayed in Fig.[T] The tt phase kink runs 
along the diameter at the direction ip — ±7r/2, and the 
phase difference saturates to and tt at the left and right 
part of the cyhnder (Fig. [TJa)). 

Associated with the static phase kink, the oscillating 




FIG. 2: (color online). Same as Fig. [T]for the (2,1) mode of 
a cylindrical mesa. 



magnetic field assumes the maximal absolute value along 
the diameter at cj) = ±7r/2, and decreases to zero at the 
left- and right- most parts of the cylinder (Fig.[T]Jd)). The 
azimuthal component of the magnetic field is always zero 
at the edge of cylinder, as imposed by the boundary con- 
dition. The electric field takes its maximal absolute value 
at the left- and right-most parts, while is reduced to zero 
along the diameter at </> = ±7r/2 (Fig. [Ijc)). 

Comparing the distributions of the tt phase kink and 
the electric and magnetic field, it is found that the (1,1) 
mode for the cylinder geometry corresponds to the (1,0) 
mode for the rectangle geometry [31]. 

The maximal value of the supercurrent takes place 
at the left- and right-most part of the cylinder for the 
present mode. The factor cos(P*(r)) from the tt phase 
kink (Fig. [TJa)) renders the supercurrent associated with 
the cavity mode Ji[Ag{r)^ positive over the cyhnder 
(Fig. mb)), which permits large bias current when the 
plasma amplitude A is enhanced at the cavity resonance. 



B. (2,1) mode for cylindrical mesa 

The (2,1) mode in the cylindrical mesa is given by 

Pi{r,t) =ujt + AJ2(^^p^ cos{2(j))sm{ujt + ^) + fiP'{r), 

(9) 

with X21 — 3.0542 the first zero of the derivative of 
J2{z). The patterns for the (2,1) mode are displayed in 
Fig. [2] There are two pairs of ±tt kinks in the azimuthal 
direction (Fig. [BJa)), which increase the frequency to 
/cji = 3.0542/a higher than the (1,1) mode where the 
TT kink running along a diameter (Fig. [Zlja)) equivalent 
to a pair of ±7r kinks in the azimuthal direction. 

The magnetic field penetrates into the cylinder along 
two directions (Fig. (Hd)), cj) = 3tt/4 and — —tt/4, and 
flows away along the two orthogonal directions, (j) = tt/A 
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FIG. 3: (color online). Same as Fig. [T]for the (0,1) mode of 
a cylindrical mesa. 




FIG. 4: (color online). Same as Fig. [T]for the (1,2) mode of 
a cylindrical mesa. 



and (p = — 37r/4, where the absolute value of magnetic 
field assumes its maximum. The electric field and the 
supercurrent become maximal at the directions of multi- 
ples oi(t> — tt's. The (2,1) mode for the cylinder geometry 
corresponds to the (1,1) mode for the rectangle geometry 

C. (0,1) mode for cylindrical mesa 

For comparison, we display in Fig. [3] the patterns for 
the (0,1) mode for cylinder geometry which has been dis- 
cussed in Ref. 31]. Since this mode given by 

Pi{r,t)^Lot + AJo{^p)^{n{ujt + ^) + fiP\r), (10) 

with Xoi = 3.8317 the first zero of the derivative of 
Jo(-z), is uniform azimuthally, the tt phase kink has 
to be compressed into the radial direction (Fig. [S^a)). 
This makes the spatial variation of the magnetic field 
steep (Fig. EJd)), and thus the wave number increases 
to k^i = 3.8317/a, which is to be excited at a voltage 
higher than the above (1,1) and (2,1) modes. The mag- 
netic field is circular in this mode (Fig. EJd)), and in 
order to match the boundary condition at the edge it 
is suppressed to zero totally. The magnetic field is also 
suppressed to zero at the center. 

D. (1,2) mode for cylindrical mesa 

The TT phase kink can occur simultaneously in both 
azimuthal and radial directions. The lowest mode is the 
(1,2) mode as shown in Fig. 2] given by 

Pi{r,t) ^ Lot + AJi{^ cos (t) sm{ujt + (p) + fiP^ir), 

(11) 



with Xi2 — 5.3314 the second zero of derivative of Ji{z). 

E. IV characteristics for cylindrical mesa 

The IV characteristics associated with a cavity mode 
denoted by the eigenfunction g{r) is well approximated 
by [33,[3il,[3i,[33] 

r {lndrgir)cosP'^ ^ 

with k the cavity frequency and Q the cross section of 
mesa perpendicular to the c axis. The prefactor of the 

0-5 I n 1 1 I , 1 1 
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FIG. 5: (color online). IV characteristics for the cylindrical 
mesa including the four lowest modes. The dimensionless 
voltage, equal to the wave number and frequency, is given 
by = xVa with x" = 1.8412, 3.0542, 3.8317 and 5.3314 
for the (1,1), (2,1), (0,1) and (1,2) mode, respectively. The 
radius of cylinder is a = 0.4 (a = 0.4At;). 
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FIG. 6: (color online). Radiation patterns for (a) (1,1) mode, 
(b) (2,1) mode, (c) (0,1) mode, and (d) (1,2) mode of cylin- 
drical mesa at the respective resonance frequency. The radius 
of the mesa is a = 0.4. 

nonlinear resonating term is 0.379, 0.345, 0.357 and 0.314 
for the (1,1), (2,1), (0,1) and (1,2) modes, respectively. 
In Fig. [5l we display the IV characteristics for a cylinder 
with the radius a = 0.4. A treatment with im pro vement 
on the peak value of current at the resonance [3l| is also 
available for the present cylindrical geometry. 



F. Radiation pattern for cylindrical mesa 

The radiation pattern of each mode can be computed 
by using the Huygens principle [33], and the result is 
displayed in Fig. [6l 



IV. ANNULAR MESA 

One of the hurdles for the present technique to excite 
EM waves of even higher frequency is the heating effect, 
since the corresponding higher dc voltage injects large 
currents into the sample resulting in severe Joule heat- 
ing. One way to overcome this effect may be to dig a 
hole in the superconductor mesa, rendering for example 
a cylindrical one to annular, which reduces the cross sec- 
tion, and thus the total current and Joule heating. The 
inner surface of an annular mesa may help leaking heat 
generated in the mesa additionally. Tailoring the shape 
of superconductor mesa however will affect the cavity 
mode, and thus the radiation frequency in a nontrivial 
way. In order to check this effect, we need to understand 
the cavity phenomenon for the annular geometry. 

A. (1,1) mode for annular mesa 

The (1,1) mode in an annular mesa is given by 



Pi{r,t) =u;t+ (AjM^p) + AnA^iI— p) I cos0sin(c^t + </?) + fiP'ip), (13) 

I 



where Ni{z) is the Bessel function of second kind, and 
ao(ai) the outer(inner) radius. The Neumann boundary 
condition should be satisfied at both the outer and inner 
surfaces, which determines the coefficient Xii for a given 
radius ratio ai/oo: 

J'iiXn-)N[ixl,) - J[{Xn)N[ix'iA) = 0, (14) 

where J[{z) and N[{z) are the first derivatives, and the 
wave number, or equivalently the frequency, is given by 
kfi = Xii/flo- The ratio of the coefficients is then deter- 
mined by 

A^/A, = -J[{xl,)/N[{xl,). (15) 

The distributions of the static phase kink, supercur- 
rent, the spatial part of the oscillating electric and mag- 
netic field are displayed in Fig. [7] for the case of ai/uo = 
1/2, with fcJi = Xii/ao = 1.3546/ao from Eq. PH). 



There are two features concerning the redistribution of 
magnetic field, comparing Fig.[7Ud) and Fig. [T]Jd). First, 
the maximum of magnetic field at the cylinder center is 
suppressed. Secondly, the magnetic field becomes along 
the radial direction in most part of the mesa when the 
center of mesa is removed, since the magnetic field should 
be normal to both outer and inner surfaces. Both of these 
two factors suppress the strength of magnetic field, and 
thus reduce the frequency. 

B. (2,1) mode for annular mesa 

The plasma term for (2,1) mode is given by 

P(r,i) = (AjJ2{—p)+AnN2{^p)) cos(20)sin(cot+(^). 

\ Go Go y 

(16) 

The wave number of this mode is determined by an equa- 
tion same to Eq. (|14p except for using J2{z) and N2{z) 
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FIG. 7: (color online). Same as Fig. [T] for the (1,1) mode 
of an annular mesa, except that the lateral coordinates are 
normalized by the outer radius and the ratio between the 
inner and outer radii is 1/2. 




FIG. 8: (color online). Same as Fig. [7]for the (2,1) mode of 
an annular mesa. 

instead of J[{z) and N[{z) as fcf^ — xii/'^o — 2.6812/ao 
for fli/oo = 1/2. The distributions of the static phase 
kink, supercurrent, the spatial part of the oscillating elec- 
tric and magnetic field are displayed in Fig. [51 which is 
to be compared with Fig. [51 Similar to (1,1) mode, the 
magnetic field becomes along the radial direction in most 
part of the mesa. 




FIG. 9: (color online). Same as Fig. Qfor the (0,1) mode of 
an annular mesa. 




FIG. 10: (color online). Same as Fig. [Tjfor the (1,2) mode of 
an annular mesa. 

supercurrent, the spatial part of the oscillating electric 
and magnetic field are displayed in Fig. [9l 

In contrast to the modes (1,1) and (2,1), the magnetic 
field is circular in the (0,1) mode. The system arranges 
the magnetic field to zero totally at the two surfaces to 
satisfy the Newmann boundary condition. Thus, the den- 
sity of magnetic flux increases when the center of cylinder 
is removed, which results in a higher frequency. 



C. (0,1) mode for annular mesa 
The plasma term for the (0,1) mode is 



D. (1,2) mode for annular mesa 

The plasma term for the (1,2) mode is 



P(r,t)= (A,Jo{^p) + AnNoi^p))sin{u;t + ip). P{r,t) = (AjJi{^p) + AnNi{^p)cos(t)smiujt + ip). 

(17) (18) 

The wave number is kg^ = Xoi/'^o — 6.3932/ao for The wave number is — Xiil^o = 6.5649/ao for 

0.1 /ao = 1/2. The distributions of the static phase kink, fli/oo = 1/2. The distributions of the static phase kink. 
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FIG. 11: (color online). Wave-number factor as a function 
of the radius ratio Oi/oo for the lowest four modes for the an- 
nular geometry, with the wave number given by fc" = 'iC jo-o- 



supercurrent, the spatial part of the oscillating electric 
and magnetic field are displayed in Fig. 1101 

There are two features in the redistribution of mag- 
netic field for this mode as displayed in Fig. [TUTd) . The 
maximum of magnetic field at the center of the cylindrical 
mesa disappears ( Fig. |D^d)), similar to (1,1) mode. On 
the other hand, the density of magnetic flux increases at 
the other two maxima of magnetic field, similar to (0,1) 
mode. 



E. Dependence of frequency on aspect ratio 

The dependence of wave number, and equivalently fre- 
quency, on the aspect ratio Oi/ao is shown in Fig.[TT] For 
the (1,1) and (2,1) modes, the wave number k°- — x^/oo 
decreases with increasing ratio ai/ao, whereas an oppo- 
site trend is seen for the (0,1) mode. The behavior of the 
(1,2) mode is a compromise of the both trends, yielding a 
minimum in wave number at ai/ao — 0.2. These behav- 
iors are well understood from the redistribution of mag- 
netic field in the corresponding cavity modes discussed 
above. 

For the (1,2) mode, the regime < flj/ao < 0.4 can 
be used to suppress the Joule heating without changing 
much the radiation frequency. Since the curves for the 
(1,1) and (2,1) modes in Fig. [TT] are quite flat, they can 
also be useful for the same purpose. The mode (0,1) and 
the mode (1,2) with ai/ao > 0.4 can be used to increase 
the radiation frequency provided the heating effect can 
be controlled. 
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FIG. 12: (color online). IV characteristics for the annular 
mesa including the four lowest modes. The dimensionless 
voltage is uj = x"'/a with = 1.3546, 2.6812, 6.3932 and 
6.5649 for the (1,1), (2,1), (0,1) and (1,2) modes respectively. 
The radii are ai — 0.2 and Oo = 0.4. 



(a) 





FIG. 13: (color online). Same as Fig. |6] for annular mesa. 
Multi reflections due to the inner surface are neglected since 
the thickness of annular mesa is very small compared with 
the EM wavelength. The radii are ai = 0.2 and Oo ~ 0.4. 



(0,1) and (1,2) mode is 0.405, 0.403, 0.396 and 0.321, 
respectively. The IV characteristics for the annular mesa 
is displayed in Fig. [T2l 



G. Radiation pattern for annular mesa 



F. IV characteristics for annular mesa 

The IV characteristics for the annular mesa is also 
given by Eq. (dH). The prefactor for the (1,1), (2,1), 



The radiation pattern of each mode is displayed in 
Fig. [13] for the annular geometry, where multi reflec- 
tions at the inner radius have been neglected since the 
mesa thickness is very small compared with the EM wave- 
length. 
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V. SUMMARY AND PERSPECTIVES 

We first address that the radius dependence of radi- 
ation frequency observed in a recent experiment using 
a cyhndrical mesa [s^l is to be understood as a clear 
indication of the right boundary condition for mesa of 
small thickness compared with the EM wavelength, with 
which the static tt kink state has been derived uniquely 
for inductively coupled Josephson junctions. Detailed 
spatial distributions of the superconductivity phase dif- 
ference with a TT phase kink are presented for various 
cavity modes of cylindrical mesa. Along with them, we 
also summarize the spatial distribution of the EM stand- 
ing waves inside the junctions which hopefully can be 
observed in experiments. 

We propose to use annular mesas to excite THz EM 
radiations. The obvious advantage of the annular geom- 
etry is the reduction of heating, since the area of sample 
is reduced which suppresses the total Joule heating and 
the inner surface of the sample may enhance heat leak- 
age additionally. The effect of removing the central part 
of a cylindrical mesa, thus rendering an annular one, on 
the redistribution of the superconductivity phase differ- 
ence, the supercurrent and the EM waves is analyzed. 
It is shown that, depending on the mode, the radiation 



frequency varies with the aspect ratio in different ways, 
and there is plenty of room for modification of the EM 
radiation by tailoring the shape of sample, while keeping 
low heating. 

Superconductor samples of annular geometry can also 
be manufactured as a waveguide resonator, where frag- 
ments of superconductor mesas are arranged in a thread 
with each satisfying the condition of small thickness com- 
pared with the wavelength of electromagnetic wave. In 
contrast with the fiber lasers developed for visible lights 
(see for example [4l|), the lasing mode propagates along 
the annular superconductor, whereas the electromagnetic 
waves within the inner surface establishes the coherent 
oscillation in thread of superconductor mesas. When the 
outer surface is shielded by wrapping a jacket made of 
appropriate materials, the wave propagating in the wave 
guide can be enhanced. 
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